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Constraining sterile neutrinos with a low energy beta-beam 
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We study the possibility to use a low energy beta-beam facility to search for sterile neutrinos by 
measuring the disappearance of electron anti-neutrinos. This channel is particularly sensitive since 
it allows to use inverse beta decay as detection reaction; thus it is free from hadronic uncertainties, 
provided the neutrino energy is below the pion production threshold. This corresponds to a choice 
of the Lorentz 7 ~ 30 for the ®He parent ion. Moreover, a disappearance measurement allows 
the constraint of sterile neutrino properties independently of any CP violating effects. A moderate 
detector size of a few 100 tons and ion production rates of ~ 2 ■ 10^^ s~^ are sufficient to constrain 
mixing angles as small as sin 26 = 10"^ at 99% confidence level. 

PACS numbers: 13.15.-(-g, 13.90.-fi, 14.60.Lm, 14.60.Pq, 14.60.St 



INTRODUCTION 

A large number of experiments have now convincingly 
demonstrated that active neutrinos, i.e. left handed neu- 
trinos which interact via W and Z exchange, can change 
their flavor. More recently, KamLAND jl[ is providing 
direct evidence for neutrino oscillations. The oscillation 
of three^ active neutrinos describes the global neutrino 
data very well, see e.g. Q. The fact that neutrinos oscil- 
late implies that at least two of the mass eigenstates have 
a non-zero mass. Most models to accommodate massive 
neutrinos introduce right handed neutrinos, i.e. states 
which do not couple to W and Z bosons. These right 
handed neutrinos can either directly provide a Dirac mass 
term or they mediate the seesaw mechanism In the 
latter case, the right handed neutrino tends to be very 
heavy with ^ 10^^ — 10^^ GeV. However, in the most 
general scenario there is a 6x6 mass matrix whose entries 
are essentially unknown. To obtain the physical neutrino 
masses at scales below the electroweak phase transition 
this mass matrix has to be diagonalized and its eigen- 
values are the neutrino masses relevant for low energy 
observations. Since the entries of the mass matrix are 
not known it can gives rise to any spectrum of eigenvalues 
and thus neutrino masses. The need to describe the oscil- 
lation of three active neutrinos implies that at least three 
eigenvalues have to be of 0{1 eV) or less. However, there 
can be — 3 additional small eigenvalues corresponding 
to light neutrino states, which due to the Z decay width 
bound would have to be sterile. Thus, we conclude that 
sterile neutrinos are theoretically well motivated by the 
observation of neutrino mass. Furthermore, we see that 
there can be one or more sterile neutrinos which are light 
enough to play a role in neutrino oscillations. Note, that 
all these considerations are entirely independent of any 
experimental claims to have seen sterile neutrinos, like. 



for example, the one by LSND [5|. 

In this work we explore how well sterile neutrinos 
can be constrained by a dedicated oscillation experi- 
ment based on a low energy beta-beam facility. For 
other physics which can be studied using a low en- 
ergy beta-beam, e.g. see [6]. The proposed experiment 
is a disappearance experiment and will be sensitive to 
Am^ — 0.5 — 50 eV^ and can probe mixing angles as 
small as sin^ 26* = 10~^ — 10~^. The oscillation proba- 
bility for our purposes is a two flavor D,, — >■ oscillation 
and the survival probability of Pe is given by 



n 1 . 2 . 2 Am L 

Fee = 1 — sm 26 sm 
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(1) 



^ We know from the invisible decay width of the Z boson 3| , that 
there are 3 active neutrinos with m < mz/2. 



The most sensitive experiments in this mass range 
looking for the disappearance of Pg have been the reac- 
tor neutrino experiments: Bugey [7'| and Chooz [8]. The 
Bugey bound is valid in the range Am^ = 0.1 — 1 eV^ 
and it can constrain sin^ 29 below 5 • 10^^. The bound 
from Chooz is sin^ 20 ~ lO"! for Am^ > 0.01 eV^ 



EXPERIMENTAL SETUP 

The concept we propose here is based on using a pure 
De beam from the beta decay of completely ionized ra- 
dioactive ions circulating inside a storage ring. The 
Lorentz boost 7 of these ions will be chosen such that the 
resulting have an energy below the pion production 
threshold. In this case, the by far most likely reaction 
is inverse beta decay on free protons. In this experiment 
we have an accurate theoretical understanding of the neu- 
trino flux, spectrum and cross section. The detector will 
be placed so close to the neutrino source that oscillation 
will happen within the detector itself. Thus the differ- 
ent parts of the detector will effectively act as near and 
far detector. This allows the cancellation of most sys- 
tematical errors in a similar fashion as in modern reactor 
neutrino experiments. 




FIG. 1: Schematic of the detector and storage ring. 



Beta-beam 

In 0, the idea of using beta decay of exotic nuclei 
to produce well defined neutrino beams was introduced. 
The basic observation is, that if a beta decaying nucleus 
is moving with a Lorentz 7^1, the isotropic neutrino 
emission will be collimated into a beam. The isotope we 
consider here is ^He, which beta decays with a half life 
of 0.81s and has an end-point energy Eq = 4.02 MeV. 
We assume that 3 • 10^^ ions are injected per second into 
the storage ring and the beta-beam complex operates for 
1.6 • 10^ s per calendar year and our experiment runs for 
a total of five years. 

We consider a magnetic field B — 5T inside the storage 
ring. Taking 7 = 30 this yields magnetic radius p = 56 m 
for the 2He^"'" ion. The useful fraction, /, of decays then 
is given by the ratio of the length of the straight section, 
S, to the overall circumference of the ring, f — „„ ^„ — . 
With 5 = 10 m, we obtain / = 2.7%. 



Event rate calculation 

Figure [T] shows a schematic of the setup we consider. 
In this scheme we will have a cylindrical detector whose 
symmetry axis is aligned with the straight section of the 
storage ring. The free parameters are (see figure [T]): the 
length of the straight section (S), the distance between 
the front end of the storage ring and the front end of the 
cylindrical detector (F), the radius (i?) and the length 
(D) of the detector. In the following, we will derive a gen- 
eral expression for the neutrino event rate as a a function 
of these free parameters. 

Neglecting the small Coulomb corrections to the 
beta-spectrum^, the lab frame neutrino beta-beam flux 
per unit length of the straight section in units of 

sr~^ MeV~^ s~^ m~^ emitted at an angle with the beam 
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FIG. 2: The un-oscillated event rate as a function of neutrino 
energy. The event rate has been calculated including all geo- 
metrical effects and with a luminosty of 3 ■ 10^^ ionss"^. The 
different lines show the result for different values of the length 
of the straight section, S, as indicated by the labels next to 
each line. 

axis is described by 

1 9 1 



47r ml f 7(1 - ^cos6») 



(Eo-E'-)E*^J{Eo-E*y 



(2) 

where rrie is the electron mass, Eq is the electron total 
end-point energy and E* is the rest frame energy of the 
emitted neutrino'^. / is the phase space factor associated 
with the beta decay of the nucleus. 7 is the Lorentz 
boost such that E* = 7iJ(l— ^cos6i), E being the neutrino 
energy in the lab frame, g = Nq/S is the number of 
useful decays per unit time per unit length of the straight 
section. 

To calculate the resulting number of events in a cylin- 
drical detector of radius R and length D aligned with the 
beam axis it is necessary to integrate over the length 5' 
of the straight section of the storage ring and the vol- 
ume of the detector. Here we assume that the beam is 
perfectly coUincar and has no transverse extension^. The 
un-oscillatcd event rate in a detector placed at a distance 
Y from the storage ring is given by 

j-S fD r-g' f-F' 

ne dx Ml de2-Ksine I 
Jo Jo Jo J E, 



dN 
dt 



dE4,^'"''\E,e)(7{E), 



(3) 



and E' 



7(l-^cosfl) ■ 

for oscillations, is 



where tanS {x,t) = y+x+t 

Note, that the baseline, relevant 
L — Y + X + £. Here, n represents the number of tar- 
get nucleons per unit detector volume, e is the detec- 
tor efficiency which is taken to be unity in our calcula- 
tion. Ejnin denotes the energy threshold for our detec- 
tion method. We work with a threshold of 25 MeV which 



We checked that these corrections are negligible for our purposes. 



^ Quantities without the '*' refer to the lab frame. 
* Note, that the beam size is of order 10~^ m, whereas all other 
length scales are of order ~ 10 m. 
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FIG. 3: This figure shows the ratio of oscillated to un- 
oscillated events as a function of the reconstructed L/E. In 
upper panel, Am^ = 10 eV^ and in lower panel Am^ = 40 
eV'^. The value of mixing term sin^ 26 = 0.1. The red (solid) 
line includes all geometrical effects and the detector resolu- 
tion, whereas the blue (dashed) line assumes a point source 
of neutrinos. 



ensures that our events are well above the backgrounds. 
cf{E) stands for the inverse beta decay reaction cross sec- 
tion which is the predominant reaction channel at 
the considered energies. Figure [2] shows the resulting un- 
oscillated event rates for S* = 10 m, F = 50 m, i? = 3.6 ni 
and D = 28.7 ni. This corresponds to a detector mass of 
1 kton. 



Optimization of the geometry 

The goal is to obtain an experimental configuration 
which has optimal sensitivity to the disappearance of 
Ue corresponding to a mass squared difference Am^ — 
1 — 10 eV^. For this optimization we fixed the detector 
mass at 1 kton, thus D is entirely determined by R or vice 
versa. We tested the values Y = {30, 50, 70, 90} m. We 
studied variations in 7 from 20 — 35. This range was cho- 
sen to stay below or close to the pion production thresh- 
old. We also changed the detector radius from 3.5—4.5 m. 



We found that within those options the configuration 
with 7 = 30, = 10 m, L = 50 m and D = 28.7 m is 
optimal. 

Using these numbers, figure [3] shows the resulting ra- 
tio of oscillated to un-oscillated event rates for two dif- 
ferent vales of Am^ as a function of the reconstructed 
L/E. The blue line assumes that the neutrinos are all 
generated at in one point in the middle of the straight 
section. The red line fully accounts for all the geometry 
effects. Clearly, for Am^ = 10 eV^ (upper panel), sev- 
eral oscillation periods can be resolved. A comparison 
between the amplitudes of those periods allows to can- 
cel systematics to a large extent (see also, figure H]). At 
Am^ = 40 eV^ (lower panel) only an average suppression 
can be observed and the sensitivity is entirely determined 
by the achievable systematic errors. 



Detector 

The detector we envisage is essentially similar to the 
MiniBooNE detector ll|, however with a cylindrical 
shape. The important background will all be beam 
based, since beam-off backgrounds will be well measured 
and cosmic ray events can be tagged with high effici ency 
as was done in similar near surface experiments 
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The beam energy has been selected to be below thresh- 
old for pion production, therefore there are few channels 
available for neutrino interactions. Only charged cur- 
rent quasi-elastic scattering on carbon and electron elas- 
tic scattering can mimic the signal of the inverse beta 
decay primary positron. These event types will experi- 
ence the same disappearance rates due to oscillations, but 
the neutrino energy reconstructed under the inverse beta 
decay hypothesis would be systematically less than the 
true neutrino energy. At these energies the cross sections 
for these background interactions are smaller by at least 
an order of magnitude, so any effect they might have on 
the measured oscillation parameters, in particular sin^ 29 
would be at most 10%, and they can be accounted for and 
corrected. 

To further reduce the impact of these beam based back- 
grounds, the detector should be optimized for the detec- 
tion of inverse beta decay. Typically this is done by tag- 
ging the primary positron with the free neutron capture 
in delayed coincidence with a mean lifetime of '~ 200 fis 
in undoped organic scintillator. However, observing the 
2.2 MeV gamma ray from neutron capture on hydrogen 
can be a significant challenge in a detector tuned to see 
events in the 50 to 200 MeV range. Additionally, the 
long delay time will put several beam bunches between 
the primary and secondary events and will therefore in- 
crease the probability of false tags. Adding gadolinium 
to the scintillator would help somewhat by increasing the 
tag energy to 8 MeV and reducing the capture time to 
^ 30 ^s. Nevertheless, even if the neutron tag was highly 
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efficient, at least some of tlie quasi-elastic events on car- 
bon will also have a correlated neutron capture tag. An- 
other possibility is to design the detector to be sensitive 
to the positron direction. We expect the elastic scatter- 
ing events to be peaked in the very forward direction, 
while the quasi-elastic carbon events will have a much 
broader angular distribution. The angular distribution 
of the hydrogen inverse beta decay events will fall some- 
where in between. Sensitivity to the angular distribution 
can be achieved by reducing the scintillation light to a 
point where Cerenkov light can be distinguished. 

To achieve an energy resolution of 10% or better over 
the 50 to 200 MeV range should be possible with a photo- 
cathode coverage of 10% if we assume approximately 
equal parts Cerenkov and scintillation light [ll|. Posi- 
tion and timing resolutions, needed to correlate events 
with beam bunches, should be achievable at the half me- 
ter and 10 ns level. With a bunch spacing of 100 ns or 
greater, this resolution would provide sufficient space be- 
tween bunches to demonstrate the rejection of non-beam 
backgrounds. At this level of detector resulution the L/E 
uncertainty is fully dominated by the unkown production 
point in the straight section of the decay ring. 



RESULTS 

Figure m shows the obtainable sensitivity in the sin^ 26- 
Am^ plane at 99% CL. Our default configuration is 
shown as red solid line in both panels. The setup pro- 
posed here improves on the existing limit for Am^ > 
0.2 eV^. In the range leV^ < Am^ < 10 eV^ the im- 
provement is one order of magnitude or better. The up- 
per panel shows how the sensitivity changes with vary- 
ing the length of the straight section. A longer straight 
section (dash-dotted line) implies a large fraction of use- 
ful decays and thus better statistics. At the same time 
the L/E resolution is reduced. As a result the sensitiv- 
ity extends to smaller mixing angles (higher statistics) 
but at smaller /S.m? (resolution). The analogous argu- 
ments holds also for a shorter straight section (dashed 
line), which yields smaller statistics but better resolu- 
tion. This improves sensitivity to Am^ > 10 eV^. Thus 
the length of the straight section can effectively be used 
to tune the experiment to the desired range of Am^ and 
it can be envisaged to have running periods with different 
straight section lengths within the same setup. The lower 
panel shows variations of the systematic error. Again our 
default setup with systematic error of 0.01 is shown as 
red, solid line. At very small value of systematic error 
of 0.001, the sensitivity (dash-dotted blue line) becomes 
essentially independent of Am^, once the first oscillation 
maximum can be observed. Note, that this very small 
value of systematic error of 0.001 probably is not attain- 
able in a real experiment. For more realistic values of 
systematic error in the range 0.005 — 0.05, the sensitivity 
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FIG. 4: Exclusion plots of sensitivity to active-sterile oscilla- 
tions in a near detector low energy beta-beam set-up. These 
are 99% CL limits (1 dof). In upper panel, we vary the 
straight section of the storage ring while in the lower panel, 
we show the performance taking various values of expected 
systematic error of the considered setup. The pink (dash- 
dot-dotted line) line shows the current, combined limit on 
1 — Pee from Bugey [Tj and Chooz @]. 



limit does not change for Am^ < 10 eV^, which is due 
to the cancellation of the normalization error between 
different oscillation maxima. This is also illustrated in 
figure O 

Figure [S] shows the obtainable sensitivity in the sin^ 2Q- 
Am^ plane at 99% CL. Our default configuration is 
shown as red solid line in both panels. The upper panel 
shows how a variation of detector mass by a factor of 10 
changes the sensitivity. Remarkably, this change is quite 
small: the detector density and aspect ratio are fixed. 
Thus with increasing detector mass, the additional de- 
tector mass will be exposed to a weaker neutrino fiux for 
purely geometrical reasons. The lower panel shows the 
change of sensitivity for a changing beam luminosity and 
here the effect is quite pronounced as every additional ion 
contributes equally to improve the sensitivity. Note, the 
our default beam luminosity of 3 • 10^^ ions s^^ is based 
on existing technology and thus may be somewhat con- 
servative. 
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FIG. 5: Exclusion plots of sensitivity to active-sterile oscil- 
lations in near detector low energy beta-beam set-up. These 
are 99% CL limits (1 dof). In upper panel, we choose different 
values of detector mass in the range 0.25 to 2 kton. Lower 
panel depicts the performance of the set-up for different val- 
ues of total injected ions. The pink (dash-dot-dotted) line 
shows the current, combined limit on 1 — Pe^ from Bugey 0] 
and Chooz [gj. 



CONCLUSIONS 

We have studied a near detector setup at a low-7 beta- 
beam facility for its ability to constrain the disappear- 
ance of electron anti-neutrinos for mass squared differ- 
ences Am^ = 1 — 10 eV^. The key point is, that for a 
suitably chosen geometry several oscillation maxima oc- 
cur within the same detector and thus a disappearance 



measurement at the sub-percent level becomes possible 
without requiring a stringent control on systematic er- 
rors. We focused on using a beam from the decay of ^He, 
which produces electron anti-neutrinos. This allows to 
use inverse beta decay as detection reaction and we can 
exploit the well defined relationship between the positron 
and neutrino energy. Thus, we have a very clean sample 
of electron anti-neutrino events. We carefully optimized 
the geometry and beam energy and found that 7 = 30 
yields the best sensitivity while still having the bulk of 
neutrinos below the pion production threshold. In order 
to have sufficient resolution in L /E we had to reduce the 
length of the straight section down to 10 m, which makes 
this setup unique. Note, that this allows our experiment 
to run parasitically in a low energy beta-beam facility 
since we use only around 3% of all ions. For a conser- 
vative beam luminosity of 3 • 10^^ ions s~^ and detector 
mass of 1 kton we obtain a sensitivity to sin^ 29 ~ 10^^ 
(99% CL) for Ato^ = 1 - 10 eV^. 
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